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Abstract 

We show that the subcritical d-dimensional nonhnear Schrodinger equation iil)t + Ai/' + IV'P'^V' = 0) 
where 1 < crd < 2, admits smooth solutions that become singular in LP for p* < p < oo, where 
P* ■= ad-i ■ Since lim(j(i^2- P* = 2, these solutions can collapse at any 2 < p < oo, and in particular 
for p = 2a + 2. 



1. Introduction 

The focusing nonlinear Schrodinger equation (NLS) 

iV't(t,x) + A^ + |V^|2-V^ = 0, V'(0,x) =Vo(x), (1) 

where x = {xi,...,X(i) G M*^ and A = Yl'j=i'^xjXj is the Laplacian, has been the subject of 
intense study, due to its role in various areas of physics, such as nonlinear optics and Bose-Einstein 
Condensates (BEC). The NLS is called subcritical, critical, and supercritical if ad < 2, ad = 2, 
and ad > 2, respectively. It is well-known that in the critical and supercritical cases, the NLS ([1]) 
possesses solutions that become singular in a finite time [Ij]. In this study we show that, contrary 
to common belief, the subcritical NLS also admits solutions that become singular in a finite time. 

Most of NLS theory has been developed for solutions that are in H^{]R.'^). In this case, the 
initial condition tpQ G H^, and the NLS solution is said to become singular at t = T^, if Tp{t) G 
for < t < Tc, and limt^Tc 11^(^)11/^1 = c«. In 1983, Weinstein proved that all solutions of the 
subcritical NLS exist globally: 

Theorem 1 ([2]). Let ip he a solution of the NLS let < ad < 2, and let £ H^. Then, 
ip exists globally in . 

Until now, this result has been interpreted as implying that the subcritical NLS does not admit 
singular solutions. In this study we show that if we do not restrict ourselves to solutions, then 
the subcritical NLS also admits singular solutions. Here, by singular we mean that there exists 
some 2 < p < oo, such that HV^Hp becomes infinite in a finite timej^ Our main result is as follows: 

Theorem 2. Let 

p*<p< oo, p*:=^—-. (2) 

ad — I 
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^In the case of solutions of the NLS, blowup of the norm implies blowup of the L'' norms for 2(7+2 < p < oo, 
see I Appendix A] 
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Then, the subcritical NLS with 1 < ad < 2 admits classical solutions that becomes singular at a 
finite time Tc in , i.e., 

< oo, < t < Tc, 



and 



lim 



Theorem [2] follows from the following Theorem: 

Theorem 3. Let p be in the range let 1 < ad < 2, let a > be a positive constant, and let 
Q{p) be the solution of 

AQ{p)-Q + ia(^^Q + pQ'^ +\Q\'^''Q = 0, < p < oo, (3a) 



subject to 
Let 

where 
and 



^ Q(0) G C, Q'(0) = 0. (3b) 
i^7''''\t^r) = -^^Qip)e-(^ (4a) 

r = |x|, L{t) = V2a(r,-i), (4b) 



1 . 1 . T, 



ds = — log— . (4c) 



L(t)' Jo L^{s) 2a ^Te-i 

Then, -0^^''^^* is an explicit solution of the subcritical NLS that becomes singular in L^ as t — > Tc . 

Remark. Although Q, hence also ip^"^^^"^, is not in H^, it is smooth, and it decays to zero as 
|x| — > oo, see Lemma [2l 
Since 

lim p* = 2+, 

then for any 2 < p < oo, there exists a singular solution of a subcritical NLS that becomes singular 
in LP. In particular, if ad is sufficiently close to 2 from below, then ■j/;^^''"*' becomes singular in 

Remark. The linear Schrodinger equation iipt + ^V' = admits the fundamental solution 
^ ~ {i-Kity/^ e^l^l^^*^*, which becomes singular in finite time in L°° Qj. Unlike ijj^"^^^"^, however, this 
solution does not become singular in L^ for any finite p. 

2. Proof of Theorem [3] 

We begin with the following result. 

Lemma 1. let ■i/;^^'"^^*' be defined as in Theorem\^ Then, ^^''^^^^ is an explicit solution of the 
NLS (P). 
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Proof. Substituting xp^^^^'"^ in the NLS (P and carrying out the differentiation proves the 
result. □ 

The result of Lemma [1] was used by Zakharov [4], and subsequently by others (see [1] and the 
references therein), in the study of singular solutions of the supercritical NLS. These solutions 
undergo a quasi self-similar collapse, in which •0^^^"^'* is the asymptotic blowup profile of the 

collapsing core of the solution. Here, in contrast, ip^^^^"^ is an explicit, "truly" self-similar solution 
of the subcritical NLS. 

We now establish the decay at infinity of all solutions of equation ([3]): 

Lemma 2. Let a > and let 1 < ad < 2. Then, for any Q{0) € C, the solution of equation ([3]) 
exists, is unique, and decays to zero as p — > oo, so that 

\Q\ = oip-^+^i"), p oo. 

Therefore, Q is in for any p* < p < oo. 

Proof. The proof is nearly identical to the proof of Johnson and Pan in the supercritical case 



see Appendix B □ 



Lemma 3. For any p* < p < oo, ip'^^^^"^ becomes singular in as t — > Tc . 



Proof. Since 



explicit / , \ 1 1 WQWp 



and WQWp < oo, the result follows. □ 
This concludes the proof of Theorem [3l 

3. The Q equation in the subcritical case 

As in [d], the far-field asymptotics of Q can be calculated using the WKB method: 
Lemma 4. Let Q{p) be a solution of equation ([3|), where 1 < ad < 2. Then, 

~ ciQi(p) + C2(52(/o), P — ^ oo, (5) 
where ci and C2 are complex numbers, and 

Qi ^ p ' , c^2 ~ e ^ ' p ' ' . 
Proof. See [Appendix C| □ 



Corollary 1. If I < ad < 2, then 

Qi G L\R'^), VQi G L\R'^), 
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and 

In addition, Qi € L^(M'^) for any 2 < p < oo, and 

Q2 £ L^iR''), < p < 00. 

ad — i 

Proof. This follows from Lemma HI □ 

In the supercritical case, a key role is played by the zero-Hamiltonian solutions of the Q equation, 
which behave as ciQi at large p [H- We now show that there are no such solutions in the subcritical 
case: 

Lemma 5. When 1 < ad < 2, there are no nontrivial solutions of the Q equation ([3]), such that 
C2 = 0, i.e., that 

Q{p) ^ ciQi{p), p — ^00. 

Proof. By negation. Assume that there is such a Q. In this case, it follows from Corollary [1] 
that Q G . Hence, ilfg^^^'^^^ is a solution of the subcritical NLS that becomes singular in , 
which is in contradiction with Theorem [TJ □ 



4. Simulations 




P P 

Figure 1: Solutions of equation ^ with d = I, a = 1.9, a = 0.5145, and with A: Q(0) = 1.2953, 
and B: Q(0) = 3. 

Figure [1] shows two numerical solutions of equation ([3]). As expected, see Lemma HI 

\Q\ ~ \cip~'^''~^/'' + C2e-'''P^ p''^-'^+^/''\ 

decreases to zero as p — > 00, while undergoing faster and faster oscillations. The "cleaner picture" 
in Figure [T]A has to do with the fact that the values of a and (5(0) were chosen so as to minimize 
the value of C2 II 



^These values of a and (3(0) were calculated using the shooting algorithm of Budd, Chen and Russel 0] for 
calculating the zero-Hamiltonian solutions in the supercritical case. 
^The value of C2 cannot be equal to zero, see Lemma [5] 
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Figure 2: Numerical solution of the subcritical NLS with the initial condition equation 'ipo{x) = Q{x) 
(solid). Dotted line is the analytic solution -j/^^''''"*. a) t = 0, L = 1; b) t = 0.4956, L = 0.7; 
c) t = 0.8163, L = 0.4; d) t = 0.9329, L = 0.2. 



In Figure [2] we solve numerically the subcritical NLS with d = \, a = 1.9, and the initial 
condition '4)o{x) = Q{x), where Q is taken from Figure [T]A.. By Lemma [U the analytic solution of 
this equation is given by ip'^^^^"^ with L = y/1 — 2at. As expected, the numerical solution agrees 
with the analytic solution, thus providing the first ever simulation of a singular solution of the 
subcritical NLS. 

In these simulations, we used a standard fourth-order finite-difference implicit scheme with dx = 
0.05 and dt = 0.001 over the spatial domain —70 < x < 70. Nevertheless, because of the slow decay 
and the ever faster oscillations as x — > oo, the agreement between the analytic and numerical 
solutions breaks down after focusing by less than 3 (see Figure [2ll). We could, of course, take an 
even larger domain with a finer mesh. In that case, the numerical solution would simply bifurcate 
from the analytic one at a higher focusing level. The point of this simulation, however, is not to 
establish numerically the existence of a singular subcritical solution (which we prove rigorously), 
but rather to illustrate the numerical difficulties in computing this solution, by showing that even 
with a relatively large domain and a fine grid, the numerical solution breaks down after focusing 
by less than 3. This suggests that numerical simulations may be useless in studying the stability 
of these solutions. 
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5. Final remarks 



Until now, it was believed that only the critical and supercritical NLS admits singular solutions. 
In this study we showed that if we do not limit ourselves to solutions, then the subcritical NLS 
also admits solutions that become singular at a finite time. This finding raises several questions, 
which are currently open. One question is whether the explicit singular solutions are stable. As 
note, this question is hard to study numerically, because of the slow decay, coupled with the ever 
faster oscillations, of the solution at infinity. Another open question is whether the subcritical NLS 
admits singular solutions that are not self-similar. The answers to these questions will determine 
whether singularity formation in the subcritical NLS will remain as an anecdote, or lead to a new 
line of research. 
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Appendix A. Blowup of ||^||p 

The NLS conserves the power (mass, norm) and the Hamiltonian, i.e., 

llV^lli = llV^olli, H{t) := - = H{0). 

a + i 

Therefore, when HV'Hhi becomes infinite, then so does ||VV'||2, hence | IV'I |2(t+2- Therefore, since 
\\ip\\2 is conserved, it follows from the interpolation inequality for norms that HV'llp also becomes 
infinite for 2a + 2 < p < oo. 
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Appendix B. Proof of Lemma [2] 

As in the proof of Jolmson and Pan in tlie supercritical case jH], let 

Q{p) = n(p)e~*'^^'/^ u{p) = ui{p) + iu2{p), 
where ui and U2 are real functions, let 

t;,(p) = p('^-i)/\(p), j = l,2, 



let 



let 



t = p^, Xj (t) = Vj (p) , yj{t) = ^ 



and let 



where 



Then, 



1 e 



Hit) = M+9i) + ^{>^-^-^]ift + /2 ) + 



1 



/3 = 1 + 



ad 



4(2o- + 2) 
1 



e = -d{d-4). 



H'{t) = ^^{hg2 - f29i) + ^ (l + t) ifi + /I) 



f3 



4(2cj + 2) 



where 



B.„(^-i,<0. 

In addition, from the Cauchy-Schwartz inequality 0], 

1 



\fi92\ < 



Since /5 > and B < 0, 



\B 



1 



H'it) < L^(|/,<^2| + |/25i|) + — 1 + - (/2 + /2). 



4t 



< 



1^1 f Vx, 



8t2 
2 



1/251 1 < 



2e\ 



1 



^f2^ 2 



8t 



Vx 



2e 



^2 . fl) 



+ /I) + ^ (5? + 5i ) ) + £o ( 1 + ^ ) (/i' + 



< 



\B 



1 



8t2 



2e\ 



* 7 



2\/At 



l + 0(-) ]H{t). 
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Since H{t) > for large t, 

H' \B\ 

Therefore, as in {5|, there exists a constant c > 0, such that 

< c(l + 0<t<oo, 



where 



Hence, 



and 

Therefore, 



LB Id 

2 a = -^L = . 

2V\ o 2 



|/,(t)|<c(l + tH), 

|^i(p)|<Cp-l/2(l + p2H)^ 

n| = 0(/)-'^+i/'^), p^oo. 



Remark. The only difference from the original proof of Johnson and Pan is that in in the 
supercritical case > 0. Therefore, we take the absolute value of i?, instead of i?, in the bounds 
for R'. 

Appendix C. Proof of Lemma [3] 

Let 

Q[p) = e-5/(^+*'^^)z(p) = e--''V4p-(rf-i)/2^(^)_ 
Therefore, the equation for Z is given by 

7in\^((^^2 1 . da-2 (d-l)(rf-3) ^ ^ 

^ (P) + I - 1 - 4^^2 + 1^1 1^ = 0- 

Since by Lemma [21 limp_^oo Q = 0, let us look for an asymptotic solution of the form 

Z = e'^^p\ w{p) ^ wq{p) + wi{p) + ■ ■ ■ . 

The equation for {wi{t)} is given by 

K + < + ..0 + K + ^ + --0^ + (^P^-i-^a^-^^^^^^ + IQl''^)=a (C.2) 

A-priori, the equation for the leading-order terms is 

< + {<? + = 0. 
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The substitution wq = cp^ shows that the order of the terms in this equation is /o"~^, /c^""^, and 
p^, respectively. Since the only consistent way to balance the leading-order terms is if n = 2, the 
equation for the leading-order terms is given by 



Therefore, 



, la I'd 2 

Wo = iyP, Wo = ±—p 



The balance of the next-order terms is given by 

, , d(j — 2 
Wq + 2wqWi — 1 — ia — = 0. 

Substituting w'q = ziziap/2 and rearranging gives, 

, i da -21 1 /zd(T-2 1\ 

ap 2a p 2p \ a 2a 2 J 

We will now show that W2 = o(l). Therefore, we obtained the two solutions 

w^'^\p) = + f-- - - logp + o(l) , 



4 V « 2 a^ 

. ^ f i —1 — d 1 \ 

w^'>{p) = -m^+ - + — — + - logp + o 1 . 
4 \a 2 a J 

Substituting Qi{p) = Q-i^pV'^ p-{d-i)/2^w(-){p) ^^^^^g ^^le result. 

In order to confirm that W2 = o(l), we note that the equation for W2 is given by 



w'i + {w[Y + 2wW2 - 4^ + \Qr = 0. (C.3) 

In the case of Qi, since IQip*^ ~ /o^^, substituting the expressions for wo and wi gives 

In the case of Q2, since |<32p'^ ~ p-2ad+2 ^ ^-2^ leading-order equation for W2 becomes 

2w'qw'2 + |Q2|^'' = 0. 

Since w'q ~ p, then w'2 ~ p-^o-d+i ^ ^-2o-d+2 ^ ^^-j.). 

Finally, we note that this this proof is rigorous, since solutions of linear ODEs always have their 
asymptotics obtained by WKB calculations, and the ODE (jC.ip for Z is "linear", since it can be 
written as 

Z"(P) + - 1 - .a'^ - ^"-'Ij^t''^ + 0(,-«/,)) z = 0. 
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